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ABSTRACT 

We present new far-infrared continuum and mill imeter- wave molecular 
observations of NGC 2024. The measurements were obtained at relatively 
high angular resolution, enabling us to describe the source energetics 
and mass distribution in greater detail than previously reported. The 
object appears to be dominated by a dense ridge of material, extended in 
the north-south direction and centered on the dark lane that is seen in 
visual photographs. Our maps of the source using the "high-density" 
molecules CS and HCN confirm this picture and allow us to describe the 
core structure and molecular abundances. The radio molecular and 
infrared observations support the idea that an important exciting star 

in NGC 2024 has yet to be identified and is centered on the dense ridge 

{/' 

about 1' south of the bright mid-infrared source IRS 2. The data 
presented here, along with other observations, allow us to describe a 
model for the source. 

Subject headings : infrared: sources — interstellar: 

molecules -- nebulae; H II regions — 
nebulae: individual (NGC 2024) — 
radio sources: general 
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I. INTRODUCTION 

NGC 2024 (Ori B, W12) is one of the best known infrared and radio 
sources in the sky. Nevertheless, because of its location close to the 
Orion Nebula region, it has received less attention than other important 
objects. Harper (1974) presented the first far-infrared map of the 
source, obtained at low angular resolution (5‘). About the same time 
Grasdalen (1974) reported detailed near- and mid-infrared observations, 
discovered a bright infrared point source (IRS 2), and suggested that it 
was the exciting star for the object. Hudson and Soifer (1976) mapped 
the object at 400 ym and compared their results with existing 
far-infrared continuum and millimeter-wave molecular data. Other 
infrared observations of note are the near- infrared spectroscopy of IRS 
1 and IRS 2 by Thompson, Thronson, and Campbell (1981), who investigated 
the excitation and structure of the source and concluded that the 
dominant exciting star had yet to be identified. More recently, Black 
and Willner (1983) discussed high-resolution near-infrared spectroscopy 
of IRS 2, emphasizing the 1 ine-of- sight CO absorption and the nature of 
the infrared source. In addition Frey et al . (1979) presented 
near-infrared continuum maps and also argued for multiple sources of 
excitation. 

Detailed high-resolution radio continuum observations of NGC 2024 
are not common. An 11 cm map has been presented by Turner ejt al . 
(1974), as well as useful single-dish measurements by Rodriguez and 
Chaisson (1978), Lobert and Goss (1978), Kriigel et al . (1982), and 
Bieging, Wilson, and Downes (1982). 

Elaborate radio molecular observations of NGC 2024 are also fairly 
rare. Bieging, Wilson, and Downes presented H^CO absorption in the 
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object and an HCN map credited to Gilmore et can be found in Hudson 
and Soifer. Recently, Snell et al . (1983) reported a detailed 
investigation of density structure in the object. In addition to this, 
there are peak antenna temperatures for a number of abundant molecules, 
but 1 ittle el se. 

In this paper we present multi-wavelength far-infrared continuum 
maps of NGC 2024 at an angular resolution of 49". We also show new J = 
2 -* 1 CS and J = 1 0 HCN maps of the source with an angular resolution 

of 80". These data are combined to discuss the excitation and structure 
of NGC 2024. In particular, we describe a fairly detailed model to 
explain the available observations of the object. This source was 
chosen as it appears we are viewing a shock/ionization front edge-on. 
Also, the relatively large optical depths in the infrared require that 
we go to longer wavelengths to study the object. Finally, since the 
exciting star has never been definitely identified, our multi-wavelength 
observations give us the opportunity to identify its location. 
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II. OBSERVATIONS AND METHODS OF ANALYSIS 
a • ) IM F ar- infrared Observations 

The NGC 2024 region was observed in the far-infrared during a 
series of flights onboard the NASA Kuiper Airborne Observatory during 
December, 1981 and January, 1983. The detector system was a 
multi-filter, seven-element array with an angular resolution of 49" for 
each detector. The array was made up of a center detector, surrounded 
by six close-packed detectors. The center- to-center separation for each 

I 

channel was 55". The reference beam spacing was 8.5 in an east-west 
direction. Ori A was used for calibration. The object was mapped both 
by scanning and with position-by-position integrations. 

The source was mapped in four filters: three narrow-band filters, 

having effective wavelengths at 40 pm, 60ym, and 100 ym, and a 
long-wavelength pass filter with a cut-on at 140 ym. The 100 ym 
observations were obtained at full -beam spacing, but the other three 
employed one-half- or one-third-beam spacing, except in the regions of 
the lowest contours. An area roughly 6' square, centered on the peak of 
the 60 ym emission, was mapped in all filters, except for the 40ym 
filter where an area roughly 4' square was mapped. As will be discussed 
below, some far-infrared mapping extended further north-south. The 
results are presented in Fig. 1. 

In addition we mapped a cross-shaped region at 200 ym with an 80" 
beam. This is presented in Fig. 2, superimposed upon a schematic 
representation of the visual appearance of NGC 2024. We do not have a 

3 

good calibration for this filter, but we estimate 1.00 = 4 x 10 Jy in 
the figure. 
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The internal uncertainty of our photometry and of the flux density 
levels in the maps is ± 0.01-0.02 (la rms) of the peak. Therefore, the 
total uncertainty of the quoted flux densities (Table 1 and below) is 
dominated by external (mainly calibration) uncertainties. From past 
experience we believe this to be about ± 20%, or less. The absolute 
positional accuracy of the maps are ± 15" (2a rms). 

Conversion of broadband fluxes to flux densities was accomplished 
in an iterative manner similar to that described in Thronson and Harper 
(1979) or Smith et al . (1982). Peak flux densities, as well as the 
total flux densities, are presented in Table 1. In a 49" beam the peak 
values are 30-40% larger than those reported by Thronson et al . (1978) 
at the same resolution. As judged by the maps in Fig. 1, this 

I 

difference is almost certainly due to the small (2.4) chopper throw used 
in the earlier observations. The total flux densities in Table 1 are 
in fair agreement with those of Harper (1974), obtained with a 5‘ beam. 
The source spectrum can be represented by an analytic function of the 
form F v “ vB v (T), where T is taken as a first approximation to the 
temperature of the dust, T^, that dominates the far-infrared emission. 
For positions removed from the peak, the same functional form was 
as sume d to characterize the spectrum, so that T^ could be found from 
observations made with only two or three filters. The temperature 
structure in the source is discussed in § III and presented in Fig. 4. 

The infrared luminosity ( 10um-200ym) of the area mapped in Fig. 1 
was 4 x 10^ L 0 at a distance of 0.5 kpc, in agreement with that found by 
Harper (1974) using a large 5' beam. This is equivalent to the total 
luminosity of an 09 zero age main sequence (ZAMS) star. As will be 
discussed below, this luminosity is significantly larger than available 
from the two well "studied embedded stars (IRS 1 and 2). 
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b.) The Radio Molecular Observations 
Observations of emission from molecules sensitive to high densities 
were made during January, 1981 using the 11 m telescope of NRAO^ at Kitt 
Peak. The 97.98 GHz J = 2 -»• 1 1 ine of CS was mapped in NGC 2024 at 1' 
spacings using a 128-channel spectrometer with 100 and 250 kHz 

The National Radio Astronomy Observatory is operated by Associated 
Universities, Inc., under contract to the National Science Foundation. 


resolution. The 88.63 GHz J = 1 -*• 0 line of HCN was similarly mapped. 
The angular resolution of the telescope at these two frequencies was 

I 

about 1.3. Absolute calibration for both objects was obtained primarily 

* ★ 

by using Ori A, for which we took T^ = 8.4 K for CS and T^ = 17.2 K for 
the J = 1^0, F = 2^"l line of HCN. Since Ori A was so close in the 
sky to NGC 2024, corrections for airmass opacity were small. The 
internal uncertainties in the observations of both molecular species are 
about 1 o = 0.3 - 0.5 K for any one position. The pointing uncertainty 
of the 11m telescope was found to be less than ± 15". The maps of the 
emission from for both molecular species are presented in Fig. 3a-b as 

'k 

the maximum T^ for the line at each position. This should be close to 
total line strength since the line width did not vary much throughout 
the source. The velocity structure, as revealed by these 
high-dipole-moment molecules, is discussed in § 1 1 1 b . 

One of the goals of this millimeter mapping program was to 
determine the HCN hyperfine ratios as a function of position in the 
source. This might then be used as a measure of the variations in 
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rnol ecu Tar abundances and column densities in the object. In NGC 2024 we 
found the average hyperfine lines to appear in the ratio 1:0. 5:0. 2 

[T*(J p = L,+ 0i> :T A< J F =1 i ^ °1 ): T A (J F = V °i )] - This is 
indistinguishible from the optically thin extreme. The discussion of 

HCN emission at the molecular peak (§ 1 1 1 b) therefore generally applies 

throughout NGC 2024. 
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III. DISCUSSION AND ANALYSIS 
a • ) The Far-Infrared Observations 
i. Morphology and Energetics 

The far- infrared observations (Fig. 1) show a resolved region of 
bright emission, extended primarily in the east-west direction at 40 pm, 
60 pm and 100 pm, and in contrast, more north-south at 160 pm. Within 
the uncertainties, the far-infrared peak is coincident with that of the 
8.4 pm emission (Grasdalen 1974) and the peak radio molecular emission 
(discussed more fully in § Illb). Note also that the higher far-infrared 
contours encompass the radio continuum peaks (Schraml and Mezger 1969, 
turner et al . 1974, Kriigel et al . 1982). We therefore suggest that the 
most likely location for an exciting star is in the general direction of 
the 60 pm maximum. This position is significantly offset from the two 
bright near-infrared sources, IRS 1 and IRS 2, which would mean that 
neither dominates the excitation of the region (see also Balick 1976; 
Frey j?t ah_ 1979; Thompson, Thronson, and Campbell 1981; and Kriigel et 
al . ). This picture is further supported by the apparently small effect 
that the two objects have upon even the lower far-infrared contours and 
dust temperature (Fig. 4 and next section). 

At longer far-infrared wavelengths (Fig. Id) and at sub-millimeter 
wavelengths (Hudson and Soifer 1976), the appearance of NGC 2024 changes 
to a more north-south elongation. This change at the various 
wavelengths is probably due to the influence of different parameters at 
the various wavelengths. At far infrared wavelengths between ^ 40 pm 
and % 100 pm, for example the emission depends very sensitively upon 
temperature, but only modestly upon emission optical depth (i.e., column 
density of warm dust). Longer wavelengths — >150 pm — are usually 
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on the Rayleigh-Jeans side of the source spectrum, and the emission is 
much less sensitive to variations in dust temperature. Column density 
variations therefore become proportionately more influential. Model 
infrared sources show this effect clearly (e.g., Leung 1976, Natta and 
Panagia 1976; Rowan -Rob in son 1979, 1980). On the basis of our 
observations then, NGC 2024 is dominated by a warm high column density 
ridge, elongated primarily in the north-south direction, toward which 
are found the exciting stars. Unfortunately our 200 ym observations 
are not complete enough to confirm this north-south extent. 

The total far- infrared luminosity of the NGC 2024 infrared source 
is close to that of an 09 ZAMS star, or an ensemble of stars of 
equivalent luminosity. The same conclusion was reached by calculating 
the ionizing flux necessary to power the radio continuum source (Krugel 
et al . 1982). This suggests that there is very little dust mixed with 
the ionized gas in the source, which is usually the case for H II 
regions of only low to moderate densities (i.e., n^ < 10^ crtT^). The 
compact or ultra-compact H II regions, on the other hand, almost 
invariably appear to have dust within the ionized zone, resulting in a 
significant reduction of the ionizing flux (e.g., Thronson and Harper 
1979; Habing and Israel 1979). The identification of a late 0/early B 
star as the ultimate power source is also in agreement with the low He/H 
ratio derived from radio recombination lines (Cesarsky 1977, Churchwell 
et al , 1978, Krugel et al . ) . 

If IRS 2 is a B0. 5 ZAMS star, as suggested by Thompson, Thronson, 
and Campbell (1981), it should make no more than about a 20% 
contribution to the total luminosity and ionizing flux of the region. 
This is consistent with our conclusions drawn from the appearance of the 
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maps in Fig. 1. Since IRS 1 was classified BO. 5 V by Johnson and 
Mendoza (1964), similar considerations apply. It is therefore not 
surprising that some of the radio continuum maps referenced above show 
modest extension in the direction of IRS 1. However, with these two 
objects making up at most roughly 40% of the total excitation, an 
additional energy source with a luminosity equivalent to about a BO Z4MS 
star remains to be found. 

An interesting feature that appears on both the 60 ym and the 
160 ym maps is a knot of emission about 4' north of the far-infrared 
peak (the 100 ym and 40 ym maps did not extend far enough north). The 
"high-density" molecular emission (§ III b) also seems to show a feature 
in this region. There is nothing prominent at this position in the very 
red photograph of Grasdalen (1974). The total far- infrared luminosity 
in a 49" beam at the position of this small feature is about 400 L , 
giving a good idea of the limiting sensitivity in our observations. 

li. Du st Temperature 

The observations at 60 ym and 160 ym were used to determine the 
dust temperature, since these two filters had the most extensive spatial 
coverage and generally the highest signal-to-noise. The results are 
presented in Fig. 4. As with most far-infrared sources, the changes in 
T d over the source are fairly gradual. This is a simple reflection of 
the very weak dependence of dust temperature upon the more significant 
source parameters: geometry, radiation field, and grain composition. 

We note that the dust temperature at, or close to, the radio molecular 
maximum is only several degrees above that deduced for the gas using 
several methods (T ^ 40-45 K; Plambeck and Williams 1979, Loren et al . 
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1981, Black and Willner 1983). This is consistent with gas heating via 

thermal coupling with the dust under conditions of at least modestly 

5 -3 

high densities (n„ > 10 cm , from the results of Goldsmith and Langer 

2 

1978; see also Snell et al . 1983). 

Only directly south of the peak does T d change rapidly. One arcmin 
south, T d ^ 35 K, and two arcmin south, £ 25 K. By contrast, T d ^ 

40 K as far as two arcmin north of the peak. East and west of the 

central position, T d changes only very little, by a few degrees. IRS 1 
and IRS 2 may contribute to the maintaining of high temperature in these 
directions. There is a hot spot \\ arcmin to the east of the flux 
density peak where T d = 58 K, but all the points surrounding this 
position are around 53 K. The most striking feature of Fig. 4 is 
therefore the very abrupt drop in dust temperature directly south of the 
60 pm maximum. Emission at the shorter wavelengths of Fig. 1, 

particularly at 40 pm, are well on the Wien side of the dust emission in 

this region. The emission is quite weak, compared to the rest of the 
source and the object therefore appears broader in the east-west 
direction. 

This temperature distribution is characteristic of exciting stars 
that are abutting a very dense molecular cloud at which we are looking 
in cross-section. In the case of NGC 2024 a dense part of the cloud 
lies just to the south of the infrared and radio continuum peaks. Rapid 
absorption of photons causes the dust temperature to drop very sharply 
over a distance of only one arcminute. Elsewhere the temperature 
gradient is quite gentle, which suggests that any extinction of photons 
in other directions is more modest. 
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This general description is supported by radio continuum 
observations (Kriigel et al . ; J. Bieging, private communication) which 
show very steep contours at the south edge of the H II region. Here the 
ionizing photons are being very suddenly absorbed. 

Our interpretation of the dust temperature map is in agreement with 
the 8.7 y m map presented by Grasdalen. Its maximum is 1' south of IRS 2 
and the emission is elongated east-west. Because mid-infrared emission 
should show the distribution of the hottest material, it is reasonable 
that the dominant exciting stars are close to the position of the 8.7 pm 
maximum. 


iii. Far-infrared Optical Depth 

With a temperature known, the 60 pm emission optical depth ( T gg) 
can be found from the flux density using F v = where fl is 

the beam solid angle. The 60 pm optical depth was thereby mapped in NGC 
2024 and is presented in Fig. 5. Much of the discussion in the 
preceding two sections can now be quantified. At far-infrared 
wavelengths NGC 2024 is very bright, which implies that the emission 
optical depth is large. Using the method described in § I la we found 
the 60 pm optical depth to be ig^ = 0.3 at the peak. Due to the 
simplifying assumptions necessary for the calculation, this value is 
probably uncertain by a factor of at least 3. Nevertheless, a high 
optical depth is also supported by observations (Lemke and Low 1972) 
that suggest optically thick mid-infrared emission and by an optical 
depth at 400 pm (Hudson and Soifer 1976) in agreement with that found 
here, if extrapolated using reasonable assumptions. 

This is fairly high for an optical depth at 60 pm, and is matched 
by similarly high values throughout the rest of the cloud. Fig. 5 is a 
contour map of each point calculated as for the peak. Although the 
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uncertainty In the absolute values for is large, the relative 
uncertainties are much lower and the figure should be a good 
representation of the distribution of far-infrared-emitting dust. It is 
not unusual to find the maximum optical depth toward the location of the 
exciting stars. However, this maximum does not fall off rapidly in 
either the north or south direction. The figure graphically shows that 
the core of NGC 2024 is dominated by an extended high-column-density 
ridge. This north-south extent is consistent with the visual appearance 
of NGC 2024 (c.f., Grasdalen, 1974, Balick 1976, Fig. 2). This 
morphology was predicted by the behavior of 160 urn emission which was 
more heavily influenced by the column density of dust than was emission 
at shorter wavelengths. 

The far-infrared optical depth can be converted to a visual 
extinction, since T g0 ^ 0.003 Ay (Gatley et al . 1977; Rowan-Robin son 
1979; Keene 1981). Therefore, Ay ^ 100 mag at the far-infrared peak, 
dropping to Ay ^ 15 mags at the outer edges of the core. It is worth 
emphasizing that the far-infrared optical depth contours do not show any 
significant effect due to the presence of IRS 1 or 2. 

On the basis of molecular observations discussed in the following 
sections, we are going to argue that the bulk of the dust is foreground 
to the ionized gas and exciting stars. In the direction of IRS 2, our 
far-infrared optical depth suggests a visual extinction close to 50 
mags. For this object Thompson, Thronson, and Campbell (1981) found only 
Ay = 12 ± 3 mags from the Brackett line series. Although there are 
significant systematic uncertainties in we consider this difference 

to be real and suggest that IRS 2 is buried only part way into the NGC 
2024 cloud as we view it from Earth. Under these conditions it appears 
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brighter in the near-infrared than would more luminous stars located at 
the far- infrared/ radio continuum peak and would thereby attract more 
attention. For example, on the basis of Fig. 5, a BO. 5 ZAMS star like 
IRS 2 (Thompson et _a_K) would appear 8 magnitude s fainter at LKJ if 
moved to the center of the H II region and behind the 
far-infrared-emitting dust. 

The various observations of NGC 2024 allow a calculation of the 

mass of the molecular cloud mapped here. A number of authors have 

derived relations between molecular column density and far-infrared 

optical depth (c.f., Evans, Blair, and Beckwith 1977; Righini-Cohen and 

Simon 1977; Sargent jet al . 1981; Keene, Hildebrand, and Whitcomb 1982; 

Thronson et a_L_ 1983; § 1 1 1 b i i ) . Because of the similarity of 

assumptions, we adopt the relation found by Thronson jit al . , t = 6 x 

10" N( CO); except for the results in Righini-Cohen and Simon, the 

other studies give similar values. Therefore, toward the center of NGC 

2024, N( 13 C0) £ 5 x 10 17 cm' 2 , or N(H 2 ) %2 x 10 23 , if L 13 C0J/LH 2 J = 3 x 

10"® (see also Black and Willner [1983J and §IIIbii). The same data may 

be used to estimate the total mass of the region mapped in the 

far-infrared. Integrating over Fig. 5 we get ^ 500 M . Averaged over 

the region mapped in HCN and CS emission (Fig. 3) this translates into 

5 -3 

densities in the range of a few times 10 cm . This mass estimate 
agrees with the virial mass derived by Snell et al. and supports the 
hypothesis of dense clumps within a more tenuous gas. The more detailed 
analysis of densities by Snell _et al . estimate about an order of 
magnitude greater values than we do here. This probably reflects the 
fact that their CS data is weighted toward higher densities. 
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b.) The Millimeter-Wave Observations 
i. Morphology and Velocity Structure 

Most of the far-infrared emission observed from regions of star 

formation usually arises from the denser, central regions of an object. 

On the basis of the preceding discussion, this appears to be the case 

for NGC 2024. For this reason we have mapped the central regions of the 

source in the emission from HCN and CS (Fig. 3). Both molecules have 

relatively high dipole moments and are therefore excited under 

conditions of high density. Emission from HCN is typically seen where 

n^ > 3 x 10^ cm” 3 (e.g., Gottlieb et aj^ 1975) and J = 2 +1 CS is 

2 5-3 

found where n u >10 cm (e.g., Linke and Goldsmith 1980, Snell et 

rL ^ — 

al . ) . Although superior in signal-to-noise, our HCN map is in good 
agreement with that of Gilmore et a]^ (1975, see Hudson and Soifer 
1976). Likewise, our CS map is in agreement with, but of superior 
quality, to that of Snell et al . 

* 

Both species show a peak in close to the location of the 
far-infrared and radio continuum peaks. The emission is generally 
elongated in the north-south direction, co-extensive with the dark bar 
in NGC 2024, the sub-millimeter emission, and the 60 pm optical depth 
distribution. The extent of emission and the fact that the molecules 
are excited by high densities also confirm the high extinction derived 
for the center of the cloud. They extend the far-infrared picture in 
important ways as well. 

Fig. 5 shows the column density of warm dust and presumably the 
associated gas. Although the molecular emission can be affected by 
temperature and radiative transfer (optical depth) effects, the 
distribution of CS or HCN emission is a good indicator of the high 
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volume density distribution. Therefore Fig. 3 suggests that the high 
c olumn density ridge in Fig. 5 is due mainly to a high volume density. 
The core of NGC 2024 is dominated by a dense, centrally-located, 
cylinder-like region. 

The preceding discussion includes at least two significant 
provisos. The first is the role of optical depth in molecular emission. 
Models (Gottlieb et aJL1975; Kwan and Scoville 1975) suggest that the 
emission from HCN is more susceptible to complicated radiative transfer 
effects than is the CS emission. Quantitative estimates of this effect 
presumably await detailed observations of higher level transitions, but 
it seems a resonable explanation for the apparent differences in 
structure between Figs. 3a and 3b. 

An additional consideration that probably affects the CS molecule 
more is that of radiative excitation. This is probably a minor 
consideration in NGC 2024 since the source is so extended and of only 
modest luminosity. This effect has been considered for NGC 2024 by 
Snel 1 et al . 

The velocity distribution in NGC 2024 shows a very interesting 


structure. Here we consider the at the line maximum of the J = 1 + 
0, F = 2 +1 line of HCN, since it was usually observed at higher 


spectral resolution (0.33 km s _i ) than was CS. As with all strong radio 


molecular anission observed from NGC 2024 (e.g., Thompson, Thronson, and 
Campbell 1981), the line center V LSR is within 1 km s _ * of +10 km s~* 
throughout the region mapped. However, just as a central ridge has been 


revealed by column and volume density tracers, so does the HCN velocity 




-18- 


There is no evidence for a north-to-south gradient greater than about 
0.3 km s" 1 over the region mapped. However, one angular resolution 
element either to the west or east of this ridge, V^ R = 10.0 ± 0.2 km 
s" 1 . Another resolution element in either direction and there is a 
further decrease to V^ SR = 9.5 ± 0.3 km s -1 . The central ridge therefore 
appears to be moving away from the surrounding parts of the cloud at 'M 
km s" 1 . 


The high-density velocity structure should be compared with that 
found for other components in NGC 2024 as this will be important 
evidence for the model described in the following section. There appear 
to be two velocity components in the object. As already noted, all 
strong molecular emission and absorption have = 9-10 km s~* 

(Thompson, Thronson, and Campbell; Bieging et al . 1982) . This velocity 
is also close to that of the strongest H I absorption (Lockhart and Goss 
1978) and C II emission (Wilson et al . 1975, Kriigel et al . 1982). 
However, the strongest ionized H and He emission is found at 5-7 km s * 
(Krugel ej: al . ) , along with weaker CO emission. Since different 
velocities may refer to different locations in space, these observations 
bear directly upon the structure of NGC 2024. 


n. Molecular Abundances and Column Densities 

Certain useful relations can be derived using the improved radio 
and infrared data presented here as well as results from the literature. 
We are first interested in estimating the column densities of the two 
molecules observed here and comparing the result with the 60 ym emission 
optical depth. The following analysis implicitly assumes that the 
far-infrared emission is co-extensive with that of the molecules. 
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Detail ed radiative transfer and statistical equilibrium models are 
required to fully analyze the HCN and CS emission (e.g., Snell et al . ) . 
Because we are satisfied with order-of-magnitude estimates, we shall 
instead adopt the assumption of thermodynamic equilibrium so that the 
molecular column density may be specified by a single excitation 
temperature. For CS we use 


1.9 x 10 12 t T (CS)A’ 
1 - exp(-7.1/T ex (CS)) 


where t 


is the optical depth of the J = 2 -*• 1 


line, T is the 
ex 


molecule's excitation temperature, andAv is the linewidth in NGC 2024 
(3.5 km s' 1 ). Under the fairly high density conditions of the core of 


NGC 2024, it is likely that T ex (CS) ^ = 45 K and that t cs -*■ 1. In 

any case, these particular values serve as upper limits for the above 

1C _? 

equation. We therefore estimate N(CS) < 2 x 10 cm ; hence, N(CS) < 6 

x 10 15 cm" 2 . Snell et__ aJL (1983) estimated the CS column density 

in NGC 2024 using a multi-level model. They calculated N(CS) = 1.6 x 
14 -2 

10 cm toward the infrared/molecular maximum. This result gives 
N(CS) = 5 x 10 14 X gQ, in agreement with the limit from our data. 

The HCN molecule is more difficult to analyze even in the 


simple manner adopted here because of significant ambiguities in 

the transfer of its radiation through a dense molecular cloud. Gottlieb 
et al . (1975) were able to fit observed HCN hyperfine line ratios with 
a farily straightforward model of optically-thin emission. Kwan and 
Scoville (1975) argue instead that significant radiation trapping can 
alter the molecule's level populations and thereby duplicate the 


appearance of optically-thin emission. However, we note that in the 
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* 

case of NGC 2024, T^(HCN) « T^, as shown in Fig. 3b. According to the 

models of Kwan and Scoville, the emission therefore has only modest 

optical depth under the high density conditions calculated for this 

object. Under these circumstances we can use the Gottlieb et al . 

calculations to estimate the HCN column density. At the far-infrared 

peak we have T*(J p = 1 ? + 0^ = 7 K, T*(J p = lg-^Oj) = 1.5 K, and T*(J p 

8 lj + 0^) = 3.1 K. On the basis of the velocity structure described 

-1 -1 

in the previous section, we shall take dv/dr = 1 km s pc . From 

13 _ 9 

Gottlieb et al . , therefore, we estimate N(HCN) = 1 x 10 cm , where we 

used n^ = 1.6 x 10 6 cm -3 from Snell e t al . . This gives N(HCN) = 3 x 

. n 13 T 2 „ -2 
10 60 cra ' 

Column densities have been calculated for a number of other species 
toward the NGC 2024. In particular, Goss et al . (1976) calculated the 
OH abundance, Lockhart and Goss (1978) found the H I absorption, Genzel 
et al . (1979) estimated both CH and H^CO column densities, Hoi 1 is and 
Churchwell (1983) studied HCO, and Black and Will ner (1983) estimated 
the amount of H^ and CO in front of IRS 2. All of these values, and 
those we estimated for HCN and CS, are collected in Table 2 where we 
present the ratio of column density to 60um optical depth. Where the 
beam size of the radio observation significantly exceeded that of the 

far- infrared data or where it referred to a position other than 

the peak corrections were made to the infrared optical depth. 

We emphasize again that the combined uncertainties of the radio and 
far-infrared data mean that Table 2 can only be an order-of-magnitude 
tabulation. It is therefore not surprising that the relative abundances 
in NGC 2024, as well as the ratios of molecular column densities to dust 
optical depth, are in general agreement with other work (e.g., the 
N(C0 )/t R relations in Thronson et al . 1983). 
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c . ) A M odel 

The observations presented here and elsewhere allow a model of NGC 
2024 to be created. We first summarize our findings. 

NGC 2024 is dominated by a high density core elongated in the 
north-south direction. The column density through the center of the 
cloud is equivalent to roughly 100 visual magnitudes of extinction. The 
peak in far-infrared emission is coincident with the maximum antenna 
temperature of CS and HCN, and close to the radio continuum maximum. 
The two bright near- infrared sources - IRS 1 and IRS 2 -probably combine 
to contribute less than about half of the total far-infrared luminosity 
of the source. An additional exciting star awaits discovery and is most 
likely located toward the radio/infrared maximum. 

In many respects our description of NGC 2024 is in agreement with 
the model suggested by Kriigel et a_K (see their Fig. 7). The H II 
region lies behind the dense molecular cloud, as demonstrated by the 
H^CO absorption studies (Bieging et al . 1982). The cause of the two 
velocity components is the acceleration of ionized gas away from a 
cavity in the cloud and in a direction somewhat toward the observer. 

Our modification to the Krugel et model is as follows. First, 
as already noted, we find no evidence for IRS 2 playing a major role in 
source excitation and we believe that this object is most likely 
embedded partway into the cloud core, as viewed from Earth (Krugel £t 
al._ were unable to satisfactorily determine the role of IRS 2 in the 
source.). Second, there seems to be very little evidence for a 
significant amount of molecular material on the backside (as viewed from 
Earth) of the H II region, except as might be necessary to confine the 
ionized gas. Visual photographs and radio continuum observations 



- 22 - 


suggest that behind the dark bar lies a diffuse - not compact - H II 
region, and our preference is therefore to minimize the amount of 
molecular material that completely surrounds it. 

On the basis of Fig. 2 we also suggest that much of the visually 
bright material seen in NGC 2024 lies forward of the bulk of the 
molecular material as viewed from earth. Our 200^&m map extends into 
the region of strong visual emission, but it shows no rapid decrease at 
these positions. This means that the dust (and gas) likewise does not 
decrease rapidly. This is difficult to quantity, however, since we have 
no dust temperature information in the outer parts of the cloud. 

We note that this picture of NGC 2024 is qualitatively similar to 
that of the M42/BN-KL region, if the latter was turned 180° away from 
our line of sight. We suggest that if we were to view NGC 2024 from its 
backside we would see a visually bright, diffuse H II region, excited by 
a centrally located BO ZAMS star. Behind this ionized zone would lie a 
very dense molecular cloud, within which is embedded a lower luminosity 
object (IRS 2), offset ~ 1 ' from the most luminous star. Except for 
differences in luminosity, and perhaps evolutionary stage, this is what 


is found in fM2. 
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TABLE 1 

Source Parameters 


Far-infrared Peak Position 

a ( 1950 )= 

5 h 39 m 13 

S ±1 S 


6(1950)= 

-1° 57' 0" 

±15" 

Flux Density 

Peak 

(49" Beam) 

TOTAL 


40pm 

8200 Jy 

33,000 Jy 


60pm 

12,000 

74,000 


100pm 

12,000 

85,000 


160pm 

5,700 

47,000 


Dust Temperature at Peak 

48±7 K 

Optical Depth at Peak (60pm) 

0.3 

Infrared Luminosity (10-200pm; 500 pc) 

4 x 10 4 L 

Mass of Molecular Cloud Core 

500 M Q 


•& 
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TABLE 2 


Col limn Densities vs 60ym Optical Depth 


X 

N(X)/t 

References for 

H 2 

*3“ 

V i 

X 

10 23 

1 

H 

7 

X 

10 21 

2 

12 C0 

3-9 

X 

10 19 

1,3 

CS 

5 

X 

10 14 

4,5 

HCN 

3 

X 

10 13 

5 

OH 

9 

X 

10 15 

6 

h 2 co 

9 

X 

10 13 

7 

CH 

2 

X 

10 14 

7 

HCO 

2-4 

X 

10 13 

8 


*Black and Willner (1983) 

^Lockhart and Goss (1978) 

^Wil son et al . (1974) 

4 Snell et aJL (1983) 

^Thi s work, § III bi i 

6 Goss et aj^ (1976); very large beam 

^Genzel _et al . (1979) 

^Hollis and Churchwell (1983) 
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FIGURE CAPTIONS 

Fig, la-d - The flux density distribution in NGC 2024 in four 
far-infrared passbands. The contours are normalized to 1.00 at the 
peak. The lowest contour level was chosen to approximately encompass 
the extent of the mapped region; the dashed lines are extrapolations. A 
typical la uncertainty is 0.01-0.02 of the peak (internal errors only). 
The small cross in the center of each map is the location of the 60 urn 
peak. The larger cross above this is the location of IRS 2. The cross 
in the west side of the figure is the location of IRS 1. The arms of 
these latter two crosses are each 15" long, the estimated 2a uncertainty 
in absolute position of the far-infrared observations. 

Fi g. 2 - The 200 y m flux density mapped in a cross-shaped pattern in 

NGC 2024. The location and beam size of the observations are 
represented as circles. The normal ization of the flux density is to 
1.00 at the peak. The extent of the visual emission from the source is 
shown as a cross-hatched distribution. The three crosses have the same 
meaning as in Fig. 1. 

Fi g. 3a-b - The corrected antenna temperature distribution for CS and 
HCN in NGC 2024. The contours are in Kelvins and the lowest contour was 
chosen to approximately encompass the extent of the mapped region; 
dashed lines are extrapolations. A typical la uncertainty is 0.3-0. 5 K 
for both molecules. The three crosses have the same meaning as in Fig. 


1 . 
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Fi g. 4 - The dust temperature distribution, determined from our 60 ym 

and 160 ym observations as described in the text. The dashed lines are 
extrapolations and the crosses have the same meaning as in Fig. 1. 

Fi g. 5 - The 60 ym emission optical depth in NGC 2024. The values were 
calculated as described in the text. The dashed lines are 
extrapolations. The crosses have the same meaning as in Fig. 1. 
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